Evolution of shear-induced melting in dusty plasma 
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The spatiotemporal development of melting is studied experimentally in a 2D dusty plasma sus- 
pension. Starting with an ordered lattice, and then suddenly applying localized shear, a pair of 
counter-propagating flow regions develop. A transition between two melting stages is observed be- 
fore a steady state is reached. Melting spreads with a front that propagates at the transverse sound 
speed. Unexpectedly, coherent longitudinal waves are excited in the flow region. 
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Applying shear can induce melting [1 5] . Experiments 
have been reported in soft materials: colloidal suspen- 
sions [(3, 0] , two-dimensional (2D) electron crystals with 
magnetic field 0, foams M, polymer glasses [l(| and 
dusty (complex) plasmas [ll| . Most of these studies were 
done with a steady application of shear. 

If shear were instead applied suddenly, one could inves- 
tigate the spatiotemporal development of shear-induced 
melting. However, there have been no such experiments 
reported in any physical system to our knowledge. 

In addition to melting, another result of sudden appli- 
cation of shear is wave excitation. When strong shear is 
applied to a 2D crystal lattice, plas tic deformation oc- 
curs, and this can cause melting [11] . When the applied 
shear is weaker and deformation is elastic, transverse 
waves (phonons) propagate through the lattice [L2| . An- 
other type of wave is longitudinal. To our knowledge, at 
least within the literature for dusty plasmas, there have 
been no reports of the excitation of coherent longitudinal 
waves due to the application of shear in either the elastic 
or plastic regimes. 

Here we seek answers to three questions. First, can 
coherent longitudinal waves be generated by applying 
shear? Second, when shear is applied suddenly, and 
melting occurs, what is the spatiotemporal development? 
Third, is there a melting front, and how does it spread? 

Physical systems that allow motion essentially on a 2D 
plane include a Wigner lattice of electrons on a liquid-He 
surface [|| , ions confined in a Penning trap 13 1 , colloidal 
suspensions 14 1 , granular fluids [lj| , vortex arrays in the 



mixed state of type-II superconductors lj|, and dusty 
plasmas levitated in a single layer 171 UM ■ 

Melting of 2D svstems has a different mechanism than 
for 3D systems [l4l: i t is sometimes termed an order- 
disorder transition [19]. Order and disorder are charac- 
terized by measures of structure. For example, the local 
six- fold bond-orientational order \<j) 6 \ @, H, H^] measures 
local order, while global order can be characterized by the 
abundance of defects 17 , [l9| , which are non- hexagonal 



Wigner-Seitz cells in Voronoi diagrams pi. Il7j|. 

Dusty plasma is partially ionized gas containing 
micron-size particles of solid matter (T3] ■ Dusty plasmas 
allow atomistic scale observation of dynamics by track- 



ing particles with video microscop y, a nd they also al- 
low laser manipulation of particles [l7j. We will exploit 
these capabilities to observe at an atomistic scale, with 
both spatial and temporal resolution, the sudden onset 
of shear-induced melting. In most previous dusty plasma 
experiments, melting was studied under steady-state con- 
ditions (by changing plasma parameters |2ll423|, va rying 
particle number 24(, and laser manipulation [ill, |25j). 
The temporal development of melting has received less 



study, and experiments that have been reported [171 126| 
relied on mechanisms other than shear-induced melting. 

Particles have an electrical charge Q and are electri- 
cally confined in a single horizontal layer where they self- 
organize with a structure analogous to a crystalline solid 
or a liquid [l7[. Coulomb repulsion is shielded with a 
screening length A n 127] . Dusty plasmas are driven- 
dissipative systems [l7J, |28[ , with frictional drag on the 
rarefied gas at a damping rate Vf [HI]. The collective os- 
cillation of particles can be characterized by the nominal 
2D dust plasma frequency uj p d p30j ]. 

For a dusty plasma suspension, shear must be applied 
differently than for most substances (e.g., colloidal sus- 
pensions 0) because the suspension does not contact a 
container. We apply shear internally within our sample, 
using laser radiation pressure ll|, |20| . 

Using the apparatus of [28] , an Argon plasma was gen- 
erated in the vacuum chamber at 15.5 mTorr, powered 
by 13.56 MHz radio-frequency voltages at 184 V peak to 
peak. The particles were polymer microspheres with a di- 
ameter 8.1 /im and a gas damping rate Vf = 2.7 s _1 [29j |. 

The particles were suspended in a single layer. Be- 
fore applying shear they self-organized in a triangular 
lattice with six- fold symmetry [17J. Particle motion was 
essentially 2D, with negligible out-of-plane displacements 
and no buckling of the particle layer. The suspension 
had a diameter w 52 mm and contained > 11 000 parti- 
cles. About 2800 particles were in the analyzed region, 
(24.7 x 20.4) mm 2 . The particle spacing was character- 
ized by a Wigner-Seitz radius [13] a = 0.25 mm. 

To apply shear stress with a sudden onset, we used 
laser radiation pressure, which applies forces internally 
within the sample [Hj]. The power of a pair of 532-nm 
laser beams was increased within 40 ms to a constant 
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level. A different constant level was chosen for each of 
our four runs, varying from 0.57 to 1.90 W per beam, as 
measured inside the chamber. The beams pointed oppo- 
sitely in the ±x directions and struck the suspension at 
a 6° downward angle. Their width and separation were 
AY; « 0.2 mm and L = 4.7 mm, respectively. To apply 
shear across the entire suspension width, the beams were 
rastered across the full suspension (at a frequency high 
enough, 200 Hz, to avoid exciting coherent waves at the 
rastering frequency (HI). 

Using a top-view camera, we recorded videos for a du- 
ration of 14 s (including 5 s before applying shear) . While 
applying shear, particles flowed out of the camera's field 
of view (FOV); they circulated around the suspension's 
perimenter and then reentered the FOV 11| . Between 
runs we turned laser manipulation off for » 20 min, so 
that the suspension had enough time to cool, solidify and 
anneal 17] . Then we calculated particle positions and 
velocities in each video frame [13] ■ 

Using the wave-spectra analysis method for particle 
random motion in an ordered lattice 3l[ , we found these 
parameters for our dusty plasma: u> p d = 72.6 s _1 , Q/e = 
-(8700±900), X D = (0.33±0.07) mm. The same method 
also yielded the transverse and longitudinal sound speeds, 
C t = (4.0 ± 0.4) mm/s and Q = (15.0 ± 1.5) mm/s, 
respectively. When coherent waves are observed experi- 
mentally, they can be identified as transverse or longitu- 
dinal waves by their speed and by the direction v x or v y 
of particle motion. 

The experiment was designed to have symmetry with 
an ignorable coordinate x. We average some quantities 
over x, as denoted by () x (averaging with cloud- in-cell 
weighting with a bin width of the lattice constant, 1.9 a). 
The kinetic energy (KE) is averaged over both x and 
y for the full analyzed region, as denoted by Q xy , We 
calculate KE = KE X + KE a , where KE X = m(v1) xy /2, 
and similarly for KEj,. Although this definition includes 
the energy of directed flow as well as random motion, we 
will report the values in temperature units (K), to allow 
comparison to other experiments. 

Results in Fig. l(a),(b) show the development of KE 
and the flow velocity after we suddenly applied shear 
stress. The KE increases dramatically, by about two or- 
ders of magnitude, Fig. 1(a). Melting occurs for laser 
powers above a threshold « 0.4 W. The flow pattern, 
Fig. 1(b), broadens until reaching a steady state at k 1 
s, when the full width at half maximum (FWHM) of the 
velocity profile is AYf =1.2 mm, which we will refer to 
as the flow regions. This width AYf is determined by a 
combination of shear stress applied by the laser, gas fric- 
tion, and shear viscosity of the suspension ll| . The KE is 
greater for velocities in the x direction than in the y direc- 
tion because shear was applied in the x direction (w 53% 
of KE X is due to directed flow in the steady state) . Scale 
lengths in the flow are ordered as AY; w a < AYf < L. 
The flow is laminar, with a Reynolds number <C 10 2 11 



In the flow region, the kinetic temperature (estimated 
from v y ) is high, with F « 20; this is well beyond the 
predicted melting point T = 160 [32j |. 

To answer our first question, we observe coherent lon- 
gitudinal waves propagating away from the flow region. 
This result might be surprising because the energy input 
was purely shear. These longitudinal waves are revealed 
in Fig. 1(c) by wavefronts with a slope corresponding to 
the longitudinal sound speed, Ci . The wavefronts emerge 
from the flow region, suggesting that the longitudinal 
waves were generated there. These wavefronts repeat, 
indicating the waves are coherent, not random. 

The coherent longitudinal waves in Fig. 1(c) have a 
period ss 4tt(jJ~J. This is the almost the same as the 
period where the dispersion relations of longitudinal and 
transverse waves cross t 33]. This observation suggests 
that the coherent longitudinal waves we observe could be 
explained by scattering from transverse waves. 

Energy is carried into the surrounding lattice by these 
coherent longitudinal waves. We will consider the speed 
of these waves below, when explaining sudden shear- 
induced melting. 

To answer our second question, we find that after sud- 
den application of shear, melting occurs in two steps. 
Examining Fig. 2(a), we can see that defects initially pro- 
liferate at a higher rate, then soon after at a slower rate, 
with a distinctive transition between. We term these 
"melting stages" 1 and 2. Afterwards, at a long time, 
there is of course a steady state. We confirmed that the 
transition is ubiquitous: it occurs for all four laser pow- 
ers we tested, and also in other suspensions with various 
values of the particle spacing. 

Our experiment provides both spatial and temporal 
resolution, which we can exploit to investigate the cause 
of the transition between the two melting stages. We 
are motivated by Fig. 1(b) to analyze time series sepa- 
rately in two spatial portions: inside and outside the flow 
regions. This yields the time series in Fig. 2(b), where 
the quantity plotted is the defect number fraction, cal- 
culated by dividing the number of defects by the number 
of particles. 

Figure 2(b) reveals that defects proliferate in differ- 
ent places at different times. In melting stage 1, they 
proliferate mainly inside the flow regions until saturat- 
ing at t = 0.16 s. Importantly, this is the same time as 
the transition in Fig. 2(a). Only afterwards do defects 
spread widely outside the flow region, as seen in the bot- 
tom curve in Fig. 2(b). Thus, the transition in the phase 
diagram arises from the different timing of when defects 
appear in the two different regions: first inside the flow 
region until defects saturate there, and then outside. 

We can quantify the physical time scales for the two 
melting stages. Candidate time scales include w"^ 1 , which 
characterizes mutual particle motion due to interparticle 
electric fields, and vj due to gas friction. Initially, when 
velocities are low, we expect friction to be small, leading 
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us to consider uj~ d as the time scale for melting stage 1. 
We find that melting stage 1 has a duration on the order 
no matter how much shear is applied, as shown 
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in Fig. 2(a). Later, as the KE increases, Fig. 2(b), dissi- 
pation by gas friction grows until it balances the energy 
input from the laser, yielding a steady state. The role of 
friction in limiting the spread of melting is suggested by 
noting that the duration of melting stage 2 is rs 1.5 vj . 

To answer the third question, we find that there is 
a distinctive melting front, and it propagates at about 
the transverse sound speed, C t . This is seen in Fig. 3, 
where the colored contours show how disorder spreads 
with time. We also draw lines with a slope correspond- 
ing to Ct] comparing these to the melting front, we see 
that they coincide. Thus, the melting front propagates at 
about Ct , not the much faster Ci . We verified this result 
using three different measures of structure. These include 
the local six- fold bond-orientational order \<j) 6 \ 0, 0, 0] 
(as shown in Fig. 3), defect number fraction, and height 
of the first peak of the local pair correlation function. 

Before this experiment, it was not obvious whether the 
melting front should spread at a rate corresponding to Ct , 
Ci, or thermal conduction. Thermal motion can be de- 
composed into incoherent waves that include both trans- 
verse and longitudinal modes. Coherent waves are also 
present, including the surprising longitudinal modes we 
found propagating away from the flow region. Thermal 
conduction due to temperature gradients 34-36| would 
also lead to a spreading of energy as time passes. 

To determine whether the flow of e nerg y is in the form 



by a distinctive transition. After defects saturated within 
narrow flow regions, they spread wider with a melting 
front that propagates at about C t . 
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of waves or thermal conduction [34H36J ]. we varied the 
laser power, which will vary the gradients. Linear waves 
propagate at the same speed, regardless of their ampli- 
tude; but energy spreads by thermal conduction more 
rapidly if the gradients are larger. We found that the 
melting front propagated at nearly Ct for all four laser 
powers we tested. Thus, we dismiss thermal conduction 
as the main mechanism for the propagation of energy 
that results in the melting front. 

Our result that the melting front propagates at about 
C t suggests that the propagation of transverse waves is 
the mechanism for transferring the energy required for 
melting in our experiment. Transverse waves, which were 
excited when shear was first applied, propagate outward 
carrying energy that can create defects. As this wave- 
front travels, it loses energy by creating disorder and by 
gas friction. Eventually, after traveling about 3 mm, the 
outward traveling wave loses so much energy that it can 
no longer generate more defects, and the steady state is 
reached. It is possible that in other systems, with greater 
applied shear, the generated coherent longitudinal waves 
might be strong enough to melt a solid lattice, so that 
the melting front would propagate at C\. 

In summary, we found that coherent longitudinal waves 
are excited in this shear-induced melting system. Apply- 
ing shear suddenly led to melting in two stages separated 
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FIG. 1: (color online), (a) Time series of the average kinetic 
energy (KE). Shear was applied suddenly by applying 0.95W 
of laser power starting at t — 0. The KE rose about two or- 
ders of magnitude until a steady state was reached at ~ 1 s. 
(b) Spatiotemporal evolution of the flow velocity {v x ) x . Two 
counter-propagating flow regions broaden to a width AYf (de- 
fined as FWHM of the velocity profile). Length scales, includ- 
ing the thickness AYj and spacing L of the laser sheets, are 
ordered as AY; ~ a < AYf < L. (c) Coherent longitudinal 
waves, generated near the flow region and propagating out- 
ward. This is revealed by the repeated pattern of wavefronts 
in this plot of (v y ) x . For comparison, the longitudinal sound 
speed Ci is indicated by the straight lines. The observed wave 
period is ~ 4jru~^ . 



FIG. 2: (color online), (a) Phase diagram, showing the pro- 
liferation of defects as the KE increases, recorded at 0.018 s 
intervals, for three laser powers. A transition between two 
stages is revealed. The quantity plotted is the defect area 
fraction; this measure of structure is the area of a Voronoi 



area, (b) Time series for KE and defects, counted separately 
inside and outside the flow regions of width AYf. The timing 
indicates that defects proliferate outside only after they sat- 
urate inside the flow regions. (The inset shows about 1/4 of 
the Voronoi diagram, for the run at 0.95 W at t = 0.16 s.) 



FIG. 3: (color online). Spatiotemporal evolution of orien- 
tational order, \(f>6\, which is defined to have a maximum of 
unity for a perfectly crystalline region HHUj. Lines arc 
drawn with the slopes Ct and Ci, starting at t = 0. The 
transverse sound speed, Ct, coincides with the melting front 
propagation. The resolution of our data along the Y axis is 
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